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1. Introduction 

In order to understand the properties of ordinary matter (baryons and mesons) it is necessary to 
understand the properties of the ground state of QCD. The best way to test the physical properties 
of a system is to vary its defining conditions in order to test its reactions. To do that in QCD we 
have to consider a system different from the vacuum but sufficiently simple in order to be able to 
study it. Some possibilities are to study extremely dense matter or matter at very high temperature 
as at the beginning of our universe. A bonus in studying the ground state of these particular systems 
is that, due to the fundamental property of asymptotic freedom, QCD simplifies a lot. In the case 
of very high temperature one expects QCD to behave as a free theory describing a non interacting 
gas of quarks and gluons. A similar conclusion would be true also at high density except that 
having to do with fermions we have to keep into account the exclusion principle. As a consequence 
a very degenerate Fermi sphere is formed and, if an arbitrary attractive interaction is present, we 
expect that the phenomenon of color superconductivity takes place ^ Qj. This is what we expect 
in QCD, since at very high density the theory can be described in terms of a single gluon exchange 
and this provides an attraction in the antisymmetric diquark state. From this we get informations 
about two asymptotic regions of the phase diagram of QCD in the variables (jJ.,T), where Li is the 
baryon chemical potential. The two regions are respectively (fx rj 0, T — > °°) and (ju — > °°, T 0). 
However we would like also to know what happens in the intermediate region and in particular we 
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would like to determine the order of the phase transitions from the hadronic phase to the quark- 
gluon plasma and to the color superconducting phase. In this context the possibility arises that 
going from the hadronic to the gluon-quark plasma phase, there is a cross-over for small chemical 
potentials and a first order for higher values of }i. In this case, the end point of the first-order line 
is called the "critical point" of QCD. The study of this intermediate region is quite complicated 
since perturbation theory cannot be applied to QCD and furthermore at finite chemical potential 
the usual lattice approach fails. In this talk I will discuss the historical path through which the idea 
of a critical point came about and some of the attempts of locating it. Therefore I will not discuss 
the physics associated to the critical point and how this point can be detected experimentally. Both 
these topics are treated by other speakers in this conference and many nice reviews about the subject 
exist in the literature, see for example [Q, f|, ^J, |], |7|] . 



2. Order parameters 

The quark gluon plasma phase can be thought of as a deconfined phase and therefore one 
would like to define an order parameter to distinguish between this and the confined phase (the 
hadronic one). Such an order parameter can be easily defined for a theory without quarks (or, 
equivalently for m q — > oo). This is the Polyakov loop ^ defined as 

L(x) =txQ.(x), Q.(x) =Pexp(i dtA (x,t)\ (2.1) 

with P = l/kT and Aq the time component of the gluon field. It turns out that the expectation value 
of the loop is asymptotically given by 

(L) « lim (2.2) 

f — >oo 

with V (r) the potential between a static quark-antiquark pair at a distance r. Therefore the confined 
and the deconfined phase are distinguished by the value of (L) 

Confined phase : (L) = 
Deconfined phase : (L) ^ (2.3) 

From a symmetry point of view, L characterizes the breaking of the center of the color group, Z(N C ), 
in the case of N c colors. From asymptotic freedom we expect that at some critical temperature T c 

(L)=0, T<T C , <L>^0, T>T C (2.4) 

For finite quark masses we expect V(r) to remain finite for r — > oo. In fact the string of the color 
flux between the two color charges is expected to break when the potential energy equals the mass 
of the lowest hadronic state, M/j. Therefore (L) does not vanish in the hadronic phase but rather 
goes exponentially to zero for M/, — > °° 

(L) w e-P M " (2.5) 
When quarks are present one can define another order parameter, the chiral condensate, (fyy), 
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Figure 1: The Polyakov looop and the chiral condensate susceptibilities in the case of two flavor QCD. This 
calculation has been made with a quark mass about four times bigger than the one needed for obtaining the 
physical pion mass [S], 

characterizing the breaking of the flavor symmetry (for instance, for three massless flavors the 
chiral symmetry would be SU (3)l <8> SU (3)r®U (l)y). In this case we expect 



{xjfy) =0 for T -> oo, (fy) / for T->0 



(2.6) 



Of course, since m q ^ 0, this order parameter never vanishes but it will have a sharp variation, or a 
crossover, close to the transition. The susceptibilities of these two order parameters 



(2.7) 
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have been evaluated on the lattice [Q] in the case of two flavors. The results are shown in Fig. 
The fi gure shows very clearly that the deconfinement and the chiral transition coincide at zero 
baryon density. 

Our final conclusion is that the phase structure is characterized by 

T < T c confined phase : (L) ps 0, (WW) ^ 
T > T c deconfined phase : (L) / 0, (iffy) « (2.8) 

Given this result, in the following we will concentrate on the chiral transition which is easier to 
deal with. 



3. First attempts to evaluate the phase diagram of QCD 

One of the first attempts to evaluate the (n,T) phase diagram of QCD was done in ref. [|io|]. 
The authors evaluated the gap equation for the chiral condensate in the approximation of one gluon- 
exchange. However the paper did not contain a discussion about the nature of the chiral transition. 



Other attempts Jl 1| , |12j ] were done using the Coulomb gauge and neglecting the retardation effects 
in the gluon propagator. This approach is simple since it is very close to a non-relativistic treatment 
and it allows to vary the static potential according to the assumptions for the gluon exchange. For 
instance, the static potential has been chosen as a 5-function, Coulomb type or confining. In all 
these cases these authors have found a second order transition in the plane (n,T). 

A completely different approach was developed in refs. [13,0]. The authors derived effective 
lagrangians for different gauge groups and a single flavor. In particular they found that for N c = 2 
the line of transition is second order, whereas for 7V C = 3 it is first order. 

The paper in ref. JT^ ] started the analysis of the problem by using a Nambu-Jona Lasinio 
(NJL) model. The case studied was Nf = 2 and N c = 3. The idea is to simulate the gluon interaction 
through an effective four-fermi coupling. Although there is no reason to expect quantitative results 



close to real QCD, one hopes that universal effects can be recovered. In ref. [15] the interaction 
lagrangian used is 

— ^sym ~\~ -^det 

Xym = \g\ [( WW? + (WiY5W) 2 + {WiYsW) 2 + (WW) 2 ] 

^det = \g2 [(WW) 2 + (WiYsW) 2 - (WiYsW) 2 ~ (WW) 2 ] (3-1) 

where J^det is the t'Hooft determinant breaking the axial symmetry U(1)a, written for Nf = 2. In 
general 

^det = \gi [det{v/(l + y 5 )w} + h.c] (3.2) 

For simplicity the authors did the choice gi = g2 = g- Therefore the model depends on 3 parame- 
ters, the mass of the quarks m, the coupling g and the cutoff A defining the model. These parameters 
can be determined at zero temperature and density by using the physical values of m n , f % and a 
reasonable value for the condensate. The results, in the plane (r,/i), are shown in Fig. ^|for two 
different choices of the parameters: I: m = 5.5 MeV, g = 5.074 x 1(T 6 MeV~ 2 , A = 631 MeV; II: 
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Figure 2: The continuous lines and the dashed lines correspond to first-order and to second-order or 
crossover transitions. The label I and II refer to different choices of the parameters, see the text. 



m = 5.0 MeV, g = 2.337 x 1(T 6 MeV~ 2 , A = 925 MeV. The corresponding values of the conden- 
sates are in the first case (xfryf) = (— 247) 3 MeV 3 and in the second one (\jf\jf) = (— 359) 3 MeV 3 . 
In the figure we see the occurrence of a critical end point where the first order transition line ends. 



Refs. [|1^, |1/7J developed an approximation scheme to QCD (today known as ladder QCD) by 
using the Cornwall, Jackiw and Tomboulis (CJT) effective action [fill]. The calculation was done 
at two loops and it is equivalent to sum up the ladder diagrams with gluon exchanged. A further 
approximation was to use an ansatz for the self-energy of the type 

I.(p,T,lx)=X(T,ii) T ^-^ (3.3) 

in order to provide an asymptotic behavior consistent with the operator product expansion. A is 
a mass scale parameter and %{T,pL) is determined by minimization of the CJT effective potential. 
The fermionic condensate is related to %(T,n) by the relation 

A 3 

(YV)t^=3 ,., X(T,h) (3.4) 
At T = /i = 0, g(T,n) is the QCD gauge coupling. The dependence of g on T and jj. was chosen 



in a way consistent with asymptotic freedom []16[]. At T = /A = the parameters were chosen to 



be A = 282 MeV, a s = 0.902. In this way (yn//), renormalized at the scale A, turns out to be 
(—197 MeV) 3 . The results of this analysis are shown in Fig. |3[ The lines denoted by Lj and 
Ljj correspond to first-order and second-order transitions (in these papers the massless quark case 
was considered). These two lines are separated by the critical end point (a tricritical point in this 
case). The dashed line y is the continuation of the points where the second derivative vanishes at 
the minimum, whereas the line 8 is the location of the points where the minima of the potential go 
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Figure 3: The continuous lines L/ and L// correspond to the first and second order transitions respectively. 
The dot is the tri critical point, designed as (ji,,T t ). The dashed lines 7 and 8 (spinodal lines) are explained 
in the text. 

from three to one (see also Fig. ||). These two lines are called spinodal lines. The regions between 
7 and Lj and Lj and 8 correspond to metastable states. We see that the qualitative results are very 



similar to the ones obtained in [15] for a completely different model. Since the tried tical point is 
where second order and first order transitions meet together, one can perform a Ginzburg-Landau 
expansion of the effective potential. This was done in [19]. By performing the expansion up to 6 th 



order in the condensate % one gets 

V(x,T^) = V(0,T,^) + a 2 (T,^)x 2 + a 4 (T,n)x A + a 6 (T,n)x 6 (3.5) 



The coefficients cij(x,T) have been evaluated in QT9Q. From this expression one easily derives the 
phase diagram in the plane (02/06,04/06)- In fact, it turns out that «6 is a positive definite quantity 
in the region around the critical point. The resulting phase diagram is illustrated in Fig. ^| (see ref. 
[20]). The phase diagram of Fig. || is obtained from this one by mapping the plane (02/06,04/06) 
into the plane (n,T), at least in the neighborhood of the tricritical point. The second order line 
corresponds to 0:2 = and 04 > 0, whereas the tricritical point is located at 02 = 04 = 0. We see 
clearly from Fig. |] the presence of the metastable regions. By using this approach it is possible to 
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Figure 4: The phase diagram arising form the Ginzburg-Landau expansion of the effective potential. The 
two dashed lines correspond to the dashed lines y and 8 of the previous figure. In particular, the line 
«2 = 0, 04 < corresponds to the line 7. The shape of the effective potential in the various regions is shown 
in the figure. 



evaluate the critical exponents around the tricritical point. In fact, when close to it we can write 



T-T r 



+ a 



, ' = 2,4 



(3.6) 



Let us introduce a quark mass term which, in this variables, is proportional to the field %, say 



V m = -hx 



(3.7) 



Then the minimum condition becomes 

h = 2a 2 % + 4a 4 x 3 + 6a 6 % 5 
and denoting by 6 either jj. or T one gets [ ]T9| ] 



(x) m 9 =o,e->e c 



(X/m q — »0, 
d(x) 

dm 1m q =Q,i 



-e c ~ m, 



1/5 



8 = 1 



7=1 



(3.8) 



(3.9) 
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where a, /3 and y are the usual critical exponents. Using these relations and the scaling relations 
for a three-dimensional system (since the finite temperature cutoff the time-like modes): 



one gets 



a = 2-3v, j3 = ^(l + Tj), 7=(2-7])v, S = ^ 



1 1 
a = -, v = -, T7=0 



(3.10) 



(3.11) 



The coefficients a, v\ and v define the behavior of the specific heat, C(6), of the correlation length, 
% (6), and of the correlation function at zero momentum, G(k —>■ 0) 



C(6) 





-a 


















-V 

















-1/2 

-1/2 



-2 



(3.12) 



(3.13) 



In 1990 when we got these results we discovered a paper by Wolff [21] showing that the two- 
dimensional Gross-Neveu (GN) model has exactly the same phase structure found by us []16{]. This 
was really interesting in view of the many similarities of the GN model with QCD. This convinced 
us that, at least qualitatively, the approximations done in our calculations did not destroy the main 
properties of QCD. Also the phase diagram with a tricritical point occurs in many physical systems. 
For instance, in the vapor-liquid transition. However, in this case the diagram usually plotted is in 
the variables (density, pressure) or (volume, pressure). What happens is that the degenerate minima 
of the first order line correspond to different densities and the line Lj splits into two different lines. 
This is shown in the (density, pressure) plane in Fig. || (see, for instance, ref. [22]). 



4. Order of the transition at zero density vs. the strange quark mass 



In 1984 Pisarski and Wilczek [23] started to investigate the order of the transition at zero 
density. In particular they investigated the dependence of the order on the number of flavors. They 
used an effective theory of QCD based on the introduction of light fields transforming as the chiral 
condensate: 

Q^YlVr (4-1) 

When close to the transition, this is a light field since the condensate, and therefore the mass term 
for <1>, vanishes at T = T c . The transformation properties of <1> under the group G = U(1)a <8> 
SU(N) L ®SU(N) R are 



<D ^ e ia U L <S>U R , U L £SU (N) l , U R G SU (N) R 
It is convenient to parameterize <I> in the form 

<Z> = <!>U, U(£SU(N) 



(4.2) 



(4.3) 
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Figure 5: The phase diagram for ladder QCD (the same as for the vapor-liquid or for the GN model) in the 
plane (l/n,p). 



In this way we separate the N 2 — 1 Goldstone fields from the condensate (p. The effective G- 



invariant lagrangian is [23] 



^=-tT{d^d^)--m 2 tx{^^)-gi (tr(&<t>)) 2 -g 2 tr(&<f>) 2 (4.4) 

To this G-invariant part a piece breaking U(1)a —> Za (N) is added 

Jzf / = c(det<£ + det<I> t ) (4.5) 

At zero temperature the symmetry breaking to SU (N)l+r is enforced through the non vanishing 
expectation value 

(<&)= 0„1 (4.6) 



In [23] the /3 -function of this effective theory has been studied and it was found that for N >3 the 
transition to the symmetric phase is first-order. The proof goes through the use of the £-expansion 
in 4 — £ dimensions and the analytic continuation to e = 1 . In this way one takes into account that, 
due to the thermal cutoff of the time-like modes, the theory is effectively three dimensional. About 
this point notice that at d = 3 the dimensions of the scalar fields are [<!>] = 1/2. As a consequence 
the sixth order term 6 is marginal and it should be included into the effective expansion when the 
coefficients of <I> 2 and <I> 4 are small, that is around the tricritical point. At zero density and close 
to the critical point T c the effective potential for the condensate (p can be taken of of the form (see 



eqs. 4.4 and 4.5) 



M 2 , X A 

T 4 > 2 +4 4 (4-7) 
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with the coefficients M 2 and A depending on the parameter of the effective lagrangian m 2 , gi, g2, c 
and N. The critical temperature can then determined as a function of the parameters by the equation 



M{T C ) = 

This result was also confirmed by lattice calculations, as shown in Table |[ 



(4.8) 



Date 



Authors 



N = 2 



N = 3 



N: 



N 



Lattice size 



1987 
1990 
1990 
1990 
1990 
1992 
1994 
1995 



Gottlieb et al. [24] 



Gottlieb et al. [25] 



Fukugita et al. [26] 



Kogut et al. [g7p 



Brown et al. [28] 



Bernard et al. [29] 



Zhu [30] 



Iwasaki et al. [31] 



crossover 
crossover 
crossover 
crossover 
crossover 
crossover 
crossover 
crossover 



I st 



V 



I st 



V 



(8, 10) 3 x4 
12 3 x 8 
12 3 x4 
12 3 x4 
16 3 x4 
12 3 x6 
(16,32) 3 x8 
12 3 (6, 18), 18 2 x 24 x (6,18) 



Table 1: A compilation of the results obtained in lattice calculations of the order of the chiral transition at 
fi = and T ^ 0. 



Following ref. [28] we may draw the phase diagram in the space of the quark masses as in Fig. 
By looking at this diagram an interesting question arises: how do we go from a second order to 




Figure 6: Solid circles is where the first order transition is seen, whereas crossover transitions correspond 
to the solid squares. The dashed circle indicates the physical point. 

a first order phase transition by varying the strange quark mass? An answer to this question was 



given in refs. [J32, 33]. The argument is the following: adding a massive quark does not change the 
effective action since the light fields are unchanged. Therefore, when close to the critical point, the 
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effective potential is still of the form given in eq. |7^. However the massive quark renormalizes 
the couplings. The resulting effect is that a variation of M 2 will change the critical temperature. 
However X could change in such a way to go through zero and become negative. If this is the case 
we know that we have to add a <p 6 term in the potential (remember that the effective theory is three 
dimensional and that such an operator is marginal). In this way we may go smoothly from a second 
order to a first order transition. 

5. Universality at non zero density 



T(GeV) 
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Figure 7: Dashed lines represent second order phase transitions, whereas the solid line corresponds to a first 
order one. The solid circle indicates the tricritical point. 

At the end of the 90's there was a big revival of the studies of the phase diagram of QCD 
prompted by the analysis made at very large density and zero temperature [0J, showing the forma- 
tion of a diquark condensate (yw) an ^ tne corresponding breaking of the color symmetry. After 
that Berges and Rajagopal [|34j] studied the coexistence of the chiral and of the diquark condensates 
in a NJL model. The phase diagram found by these authors is shown in Fig. ^ and it shows the 
presence of the tricritical point. The authors justified the presence of the tricritical point by using 



an argument very close to the one used in refs. Q32| , p3fl in the case of the strange quark that we 



have discussed in the previous section. The idea is that at zero quark mass the theory belongs to the 
0(4) universality class (Ising Z2 for m q 7^ 0) and this is not changed at finite density as shown in 
ref. fl35l]. However the renormalization of the coefficients in the effective action due to the presence 
of the chemical potential might change the coefficient of the quartic term forcing the introduction 
of a sixth order term in the potential. Again, the 6 th order term gives rise to a tricritical point in the 
phase diagram. Also, recalling that the 6 operator is marginal, one expects, at most, logarithmic 
corrections to the critical exponents evaluated before. 

After the previous paper many authors reconsidered the problem of QCD at finite density and 
temperature using many different approaches. We will give here a brief list of papers delaing with 
the problem. In 1998 Halasz et al. [t36l] considered a random matrix model, for the two-flavor case, 
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Figure 8: We show the first order lines as curves at constant quark mass, m. At m = the second order line 
is shown. For m ^ the first order line ends into a critical end point. 



In ref. [ J37[ ] the chiral phase transition has been examined both in a linear a -model and in a 
NJL model. The two cases are illustrated in Fig. |9[ Once again the results agree very well with the 
universality arguments. 
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Figure 9: The two panels illustrate the phase diagram for the ff-model (left panel) and for the NJL model 
(right panel). The middle curves are the critical lines whereas the outer lines are the spinodal lines. 



Another different approach was considered in ref. [38]. The calculation was done within the 



context of the statistical boostrap principle, and again it agrees with the universality hypothesis, see 



Fig. 10. It was also found that the critical chemical potential is non zero for a large range of values 



of the bag constant B (B 1 / 4 < 282 MeV). 

As a last example we report a recent calculation made by Hatta and Ikeda [ 39 ] using again 
ladder QCD with the help of the CJT potential as in refs. Jl6| , O ). This calculation takes also into 
account quark masses and it shows the dependence of the critical end point with this variable, as 



shown in Fig. 1 1 



A comparison of the locations of the critical end point evaluated in different models can be 
found in a recent review by Stephanov This comparison is particularly interesting since it 
shows that, although different models agree qualitatively well, from a quantitative point of view 
they are quite different. For instance, at the critical end point, the critical value of the chemical 



13 



QCD critical point: 



i i i | i — i i | i i i — | i i ■ i — n — r-p — I i I ■ I — I I ■ i i I — I — i — I — n— 




iu: 20u ;uu luu suu slju , l uV BVu -uu iuuu 
Burvuji CtiLjnL^tll Potential. LI. (MeV) 



Figure 10: The (ju, T) phase diagram from the statistical boostrap model. B is the value of the bag constant. 

120, ^ , , , , , . , , , 




u [MeV] 

Figure 11: Quark masses are evaluated at the momentum scale 1 GeV. The solid and dotted lines correspond 
to first order and second order phase transitions respectively. The solid and open circles denotes the tricritical 
and the critical end points. 

potential varies between roughly 300 MeV up to about 1000 MeV, whereas the critical value of the 
temperature goes between 40 and 170 MeV. Clearly in order to have a quantitative improvement 
one would need a first principles calculation. 

6. Lattice calculations 

As noticed in the previous section one would really need to have the possibility of testing on the 
lattice the phase diagram of QCD. However the usual sampling method, based on a positive definite 
measure in the euclidean path integral, does not work in presence of a real chemical potential, since 
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the fermionic determinant is then complex. In fact, let us define euclidean variables through the 
following substitutions: 

xq->-v&, x ! '^4> 7o^7e, Y^-iyk (6-1) 
The euclidean Dirac operator in the presence of a chemical potential is 

D(n) = 7%D% + nrl, D^dg+iA^ (6.2) 

At jU = the eigenvalues of D(pL) are pure imaginary and also, if \X) is an eigenvector of D(0), 
then 75(A) belongs to the eigenvalue — A5, as it follows from 

D(O) 1 " = -D(0), 7 5 D(0)7 5 = -D(0) (6.3) 

Therefore 

det[D(0)] = n( A )(~ A ) > ^ 6 - 4 ) 

A 

At jit 7^ this argument does not hold and we lack the positivity property. However, if one considers 
the chemical potential associated to the isospin, since this is related to the conserved current T3, the 
positivity can be proved by using Ti in conjunction with the hermitian conjugation. 

Recently there have been numerous different attempts to improve the lattice calculations at 



Rewighting method (see, for instance [40, 41, B3]). 



• Taylor expansion for small /I (see, for instance []43|], Q44]], [}45j]) 

• Imaginary chemical potential (see, for example ®, @ @, @]) 
6.1 Re weigh ting 



The reweighting technique (for a review see ref. Q51|])is based on the following identity for the 
partition function 



J1 det[D(0)] \<tet[D(0)]/„^, 

Since the integration measure is taken at /I = it is positive definite. However in the numer- 
ical calculation problems arise. The ratio of the two determinants oscillates and there are large 
cancellations. Also, since the reweighting corresponds to the ratio of two partition functions with 
different actions, it decays exponentially according to the difference of the free energies, AF. This 
is proportional to the volume and therefore the statistics required for a given accuracy increases 
with the volume. An improvement of this technique is the so called "multiparameter reweighting" 
which is a generalization of the previous method pO]]. The idea is to reweight also in the lattice 
gauge coupling, writing 

z= /i^M\ (6 . 6) 

The second reweighting parameter can be used to allow the statistical ensemble to fluctuate between 
the phases and to avoid that the ensemble goes away from criticality. One of the most recent 
calculations using this technique was done in ref. [[4^] and the result is shown in Fig. 12 
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Mb (MeV) 

Figure 12: The most recent determination with the reweighting procedure is given in ref. [^2|. The errors 
are due to the reweighting procedure and on the error of the scale determination at T = 0. The values of the 
critical temperature and chemical potential are T = 162 ± 2 MeV , pi = 360 ± 40 MeV. 

6.2 Taylor expansion 




Figure 13: A comparison of the critical line in different simulations, FP Q, FK All @, EL @. 
The box is the position of the critical end point evaluated by Fodor and Katz in [pl||. 

This method makes use of the multiparameter reweighting and at the same time of a Taylor 
expansion in the chemical potential. Although this method is not useful for determining the critical 
end point, it is of interest in the heavy ion physics where values of /J. of a few ten MeV are important. 
What one does is simply to expand in a Taylor series of pL /T the reweighting factor. In particular 
this method can be used to evaluate the behavior of the critical line at small '. This method has 



been used in ref. [43] for the two-flavor case. The results are shown in Fig. ttl 
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6.3 Imaginary chemical potential 



If jJ. is pure imaginary the fermion determinant is positive and numerical simulations can be 
done easily as for the case pL = 0. Using the fact that the observables are analytic functions of ji 
except that on the critical line, one computes expectation values at imaginary jj. and then one fits 



them by a truncated Taylor expansion [46, 47]. Some of these results are given in Fig. 13 



7. Isospin chemical potential 



To end this review we will report also some result obtained in presence of an isospin chemical 
potential, /I/. The case of ^ is, in principle, interesting for the heavy ion physics. It is also 
interesting since for jx = and jX[ 7^ the fermionic determinant is positive [|^, |53|]. The problem 
(jit and jXj 7^ 0) has been studied using effective lagrangians [54, 55], random matrices [|5"6|], NJL 



model \ pl[ |58[ ] and ladder-QCD [p9[]. The most interesting effect in these studies appears to be the 
splitting of the critical line and of the critical end point. In fact, this effect could bring down the 



end critical point to a region more accessible to heavy ion experiments. We show in Fig. 14 the 



result in the NJL model [p7[]. The result is qualitatively compatible with an analogous calculation 
made in the ladder-QCD model [p9|]. However a too strong mixing of the up and down flavors 



could destroy this interesting result, as shown in ref. [50] 



T (MeV) 



= 



T (MeV) 



fi T = 30 MeV 



(uu), (dd) 



(uu), (dd) 



Ms (MeV) 




100 200 300 400 
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Figure 14: The phase diagram for the NJL model as studied in ref. 
\ii = 0, whereas in the right panel jXi — 30 MeV. 



7J. The left panel shows the case 



8. Conclusions 

As we have shown there has been a lot of activities and results in our understanding of the 
phase diagram of QCD. However most of the progress is still at a very qualitative level. We would 
like to be able to locate the critical point with a good accuracy, but for that, a breakthrough is real 
necessary. This could come by devising some clever technique in lattice calculations, or, may be, a 
new analytical method. 
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Figure 15: The phase diagram for the ladder-QCD model as studied in ref. p^]. 
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